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Abstract: It is well known that thermally pulsing Asymptotic Giant Branch stars with low mass play 
a relevant role in the chemical evolution. They have synthesized about 30% of the galactic carbon 
and provide an important contribution to the nucleosynthesis of heavy elements (A>80). The relevant 
nucleosynthesis site is the He-rich intermediate zone (less than 10~^ Mq), where a(2a, 7)^'^C reactions 
and slow neutron captures on seed nuclei (essentially iron) take place. A key ingredient is the interplay 
between nuclear processes and convective mixing. It is the partial overlap of internal and external 
convective zones that allows the dredge-up of the material enriched in C and heavy elements. We 
review the progresses made in the last 50 years in the comprehension of the s process in AGB stars, 
with special attention to the identification of the main neutron sources and to the particular physical 
conditions allowing this important nucleosynthesis. 
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1 Introduction 

For a long time our understanding of the s-process nu- 
cleosynthesis and stellar evolution advanced in paral- 
lel, but only recently they merged in a unique theory. 
The two stories started more or less at the same time, 
namely in the second half of the '50, with the first at- 
tempts to identify the main nuclear processes respon- 
sible for the synthesis of the elements beyond iron and 
with the first comprehensive picture of the physical 
properties of post-main-sequence evolutionary phases. 
It was during the '70 that thermally pulsing Asymp- 
totic Giant Branch (AGB) stars were identified as a 
promising astrophysical site for the s process, but only 
in the last decade of the XX century it was realized 
how the so called main and strong components of the 
cosmic s-process (including elements between Sr and 
Bi) can be produced in the He-rich intershell of low 
mass AGB stars. In this paper, we critically analyze 
the main steps of these two stories up to the most re- 
cent attempts to reach a unified theoretical scenario. 

2 The beginning: 1950-1960 

In 1951, Bidelman & Keenan firstly proposed a "par- 
allel" spectral sequence of giant stars, including S, R, 
N types as well as a new possible Ba II class. Stars 
belonging to these spectral classes show anomalously 
large abundances of certain heavy elements. One year 
later, Merrill (1952) announced the discovery of tech- 
netium in S stars. It was a very important discovery, 
because Tc has no stable isotopes and, for this rea- 
son, its detection suggested that the synthesis of this 
element is an ongoing process in these stars. 

In these same years, Sandage & Schwarzschild (1952) 
carried on the first models of Red Giant Stars and 
Salpeter (1952) proved that the a-\- a ?i*Be followed 
by **Be-fa — >^^C processes provide the way for He 



burning in giant stars. The theory of stellar evolution 
beyond the main sequence was born. 

The Cameron (1955) paper entitled "Origin of anoma- 
lous abundances of the elements in Giant Stars" con- 
tains the first suggestion that neutron captures by iron 
seed nuclei are a fundamental nucleosynthesis process. 
He identified the "C(a, nf^O as a promising neutron 
source in giant stars. A naive evolutionary scenario 
was depicted: as a by-product of shell H burning via 
the CNO cycle, some ^"^C is left in the hot He-rich core 
of giant stars, where the ^^C{a, n)^®0 reaction and the 
consequent s-process nucleosynthesis take place. How- 
ever, Fowler, Burbidge and Burbidge (1955) argued 
that too few neutrons would be available in that case 
(namely 1.4 neutrons per seed), because of i) the rela- 
tively low amount of ^^C in the ashes of the CNO cy- 
cle and ii) the strong neutron poisoning of ^''N, whose 
neutron capture cross section is particularly large. A 
further problem concerned how these heavy elements, 
once produced in the stellar core, become observable 
at the stellar surface. In order to overcome these prob- 
lems, Cameron (1957) for the first time introduced the 
need for an interplay between mixing and nucleosyn- 
thesis. In particular, he claimed that mixing of a few 
protons could take place between the H-rich envelope 
and the He-rich core of low mass stars during the so 
called core He fiash. This mixing would allow the pro- 
duction of fresh ^^C (via the i2C(p, 7)"N(/3-!/)"C 
chain) in the core and, at the same time, the dredge- 
up of heavy elements. Actually, such a mixing occurs 
only if the entropy barrier set up by the H burning 
shell is switched off, as happens in extremely metal 
poor stars (e.g. HoUowell, Iben & Fujimoto 1990, Pi- 
cardi et al. 2004). In the same paper, Cameron also 
noted that since the neutron capture by ^^N mainly 
proceeds through the p channel, rather than through 
the 7 channel, the neutron poisoning effect of the '^^N is 
partially reduced to a proton recycling. These protons. 
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indeed, induce an additional production of ^'^C. Figure 

2 of the Cameron paper is particularly enlightening. It 
shows how a substantial overproduction (10^ times the 
initial abundance) of light-s elements (Is), like Sr, can 
be obtained with a relatively low number of neutrons 
per seed (about 10), while the same overproduction of 
heavy s elements (hs), like Ba, and Pb require about 15 
and 20 neutrons per seed, respectively. As a matter of 
fact, only recently have we understood that the strong 
component of the s process most likely comes from low 
metallicity AGB stars, where the lower amount of iron 
favors the production of lead (e.g. Busso, Gallino & 
Wasserburg, 1999). 

Another important indication can be found in Cameron'i 
seminal paper (1957b). When discussing the spectra 
of R Andromedae, an S star, he explained how the de- 
tection of strong niobium lines, unlike the detection of 
Tc lines, is evidence that neutron capture nucleosyn- 
thesis is not an ongoing process in this star: indeed, 
®^Nb, the only stable isotope of niobium, is destroyed 
by neutron capture during the AGB, but, later on, it is 
produced by the delayed ^^Zr decay (half- life 1.6 x 10^ 
yr). Hence, an appreciable fraction of this half- life has 
passed since neutron production ceased in the material 
that has been mixed to the surface of R Andromedae. 
Later on, Cameron (1960), proposed ^^Ne{a,n)^^M.g 
as a possible alternative neutron source. 

3 The golden age of the clas- 
sical analysis: 1961-1980 

A fundamental step ahead was taken by Clayton (1961) 
(see also Seeger, Fowler and Clayton 1965), who demon- 
strated that a superposition of several neutron expo- 
sures is required to reproduce the experimentally ob- 
served abundance distribution of the s-process isotopes 
in the solar system. As an example, Seeger et al. 
proposed an exponential distribution of neutron ex- 
posures. This simple mathematical law was widely ac- 
cepted for more than 20 years as the paradigm of the 
so called s-process classical analysis. 

Meanwhile, the availability of new computer facil- 
ities in the USA, Japan and Europe allowed the de- 
velopment of accurate models of the advanced phase 
of stellar evolution, in particulax the AGB phase. So, 
Weigert (1966) and Schwarzschild & Harm (1967) dis- 
covered the occurrence of thermal pulses in the AGB 
phase. Powered by the periodic ignition of He shell 
flashes, a series of convective zones mix the products 
of the 3a reactions (essentially carbon) through the 
He-rich layer. Ulrich (1973) noted that, if a potential 
neutron source is activated during a thermal pulse, the 
partial overlap of successive convective shells naturally 
leads to an exponential distribution of neutron expo- 
sures, precisely as suggested by Clayton. 

A simple mathematical law appeared appealing for 
the classical analysis and, in addition, AGB models 
provided a natural support to this choice. Then, by 
applying Occam's razor, the paradigm of the s process 
was definitely established. The astrophysical site of 
the s process was identified in the convective He-rich 
zone of thermally pulsing AGB stars. The resulting 



thermal energy is, in this case, about 25 KeV, corre- 
sponding to the typical maximum temperature at the 
bottom of the convective zone generated by a thermal 
pulse (about 3 x 10^ K). Hence, nuclear physics inves- 
tigations were concentrated on this energy range (see 
Bao and Kappeler 1987). 

4 The search of the main neu- 
tron source, successes and 
failures: 1980-1993 

The last piece of the puzzle, the neutron source(s), was 
still lacking. Meanwhile, the classical analysis grew 
and significant progress was made in modeling AGB 

stars (Iben 1975, Sugimoto & Nomoto 1975, Truran & 
Iben 1977, Iben 1981). In particular, it was realized 
that when the thermal pulse starts, ^^Ne is synthesized 
from the ^'^N left behind by the H- burning during 
the interpulse period, after two a captures and one /? 
decay. At its maximum size, the convective zone in the 
He-rich intershell extends from the position of highest 
He-burning rate to just below the external border of 
the H-exhaustcd core, but no further. Then, enough 
neutrons could be released by the ^■^Ne(a, n)'^^Mg, but 
only if T > 3.5 x 10* K. This condition is achieved 
in massive AGB only (M > 5 M©). Following the 
disappearance of the convective shell, the base of the 
convective envelope extends down to the outer portion 
of the region previously contained into the convective 
shell and allows the mixing of the material enriched 
in C and heavy elements within the envelope (the so 
called Third dredge-up or TDU). 

However, the fact that the ^^Ne(a, n)'^^Mg is very 
marginally active in low mass AGB stars of nearly solar 
metallicity, because the maximum temperature in the 
He-rich zone is too low, is at odds with the evidence 
that the bulk of the s-enriched AGB stars in the disk 
of the Milky Way are low mass stars (Iben 1981). The 
prediction of extant stellar models is that less than 
1 neutron per '"'^Fe is produced by the ^^Nc neutron 
source! The final word against the ^^Ne neutron source 
was pronounced by Busso et al. (1988). On the base of 
stellar models of massive AGB stars, as computed by 
means of the FRANEC code in the version described 
by Chieffi and Straniero (1989) they concluded that 
the reaction ^^Nc(a, n)^^Mg is unable to give rise to a 
solar distribution of s-nuclei. 

An alternative neutron source capable of providing 
a suitable neutron flux even in low mass AGB stars, 
was actually searched for since the beginning of the 
1980s. Sackmann (1980) suggested the possible occur- 
rence of proton ingestion into the He convective shell 
at its maximum size, but Iben (1982) showed that no 
mixing of protons occurs during the thermal pulse be- 
cause at that time the H-burning shell is very efficient, 
producing a substantial entropy barrier. On the other 
hand, Iben & Renzini (1982) noted that, as a conse- 
quence of the huge expansion of the envelope caused 
by the He-shell flash, a few hundred years after the 
disappearance of the He-rich convective zone the H- 
burning shell dies down, a condition more favorable for 
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the downward proton diffusion. They proposed that, 
at that epoch, a semiconvective layer may form, driven 
by an (assumed) peak of the carbon opacity axound 
T = 10 K and that the resulting (partial) mixing may 
induce a certain proton leak from the envelope down 
to the top of the He-rich intershell. Later on, when 
the temperature at the top of the He-rich intershell 
rises, a ^^C pocket may form. In the Iben and Ren- 
zini original scenario the s process occurred when the 
^^C pocket is engulfed into the convective zone gener- 
ated by the subsequent thermal pulse and it is burned 
at relatively high temperature through ^^C(a, n)^^0 
(see also HoUowell and Iben 1989). In this case, up 
to 26 neutrons per ''^Fc arc available at solar mctallic- 
ity and the neutron density is about lO^'-lO^'^ neutrons 
per cm^. The first complete s-process nucleosynthe- 
sis calculation based on stellar inputs was obtained by 
Gallino et al. (1988). By adopting the HoUowell and 
Iben model, they showed that the ^^C(a, 'n,)^''0 reac- 
tion provides, at variance with the '^■^Nc(o, ri)'^"Mg, a 
neutron source capable of reproducing the solar system 
distribution of s-only isotopes with ^ > 80. 

In spite of these successes, the crisis was, once 
again, knocking at the door. The advent of high- 
resolution spectroscopy era allowed powerful tools to 
probe the physical conditions at the s-proccss site (Smith 
& Lambert 1986, Malanoy 1987, Malaney & Lambert 
1988, Busso et al. 1995, Lambert et al. 1995, Abia 
et al. 2001). These studies used important branch- 
ing ratios sensitive to the physical condition at the 
s-process site, like the temperature or the neutron den- 
sity (see the Kappelor et al. contribution to this vol- 
ume). In particular, the Rb/Sr ratio, which is influ- 
enced by the branching at *^Kr, is sensitive to the 
neutron density. These studies showed that the low 
Rb abundance, compared to Sr, as found in many S, 
C and Ba stars, implies a quite low neutron density 
(less than 10* neutrons per cm^). This is incompati- 
ble with any s-process scenario in which the neutron 
irradiation takes place in the convective zone gener- 
ated by the thermal pulse. Similarly, the low value of 
the solar system isotopic ratio ^^Zr/^^Zr indicates that 
the main s-process nucleosynthesis should occur at low 
neutron density. 

Also the stellar theory was reaching a (icad-eud. 
Bazan & Lattanzio (1993), investigating the energetic 
feedback occurring when the '^^C pocket is engulfed 
into the convective shell and the s-process takes place, 
found likely substantial alterations of the thermal pulse 
characteristics and those of the related nucleosynthe- 
sis. In particular, they found a high neutron density 
(10^" cm^^), implying an overproduction of both Rb 
and ^^Zr, in clear contrast with the observation. 

5 Toward a unified approach 
to stellar evolution and nu- 
cleosynthesis in AGB stars 

The solution of the stellar s-proccss enigma was found, 
perhaps by chance, by Stranicro et al. (1995). The aim 
was to check, by means of the FRANEC code, the dra- 



matic result obtained by Bazan and Lattanzio. Follow- 
ing the original suggestion by Iben and Renzini (1982), 
Straniero et al. (1995) assumed that after the 14th 
dredge-up episode of a 3 Mq model of solar metallicity, 
a small number of protons (10~® M©) reached the top 
of the He-rich zone. They then followed the evolution 
of these protons through the whole interpulse period 
and up to the development of the subsequent ther- 
mal pulse. As expected, at H reignition a ^^C pocket 
rapidly formed, but during the interpulse large enough 
temperature were developed for the ^^C(q, n)^®0 re- 
action to be activated well before the end of the inter- 
pulse. In practice, the ^^C burning and the s process 
occur in radiative conditions and before the onset of 
the subsequent thermal pulse. The neutron density is 
only 10^-10^ cm~^ and the typical thermal energy is 
just 8 KeV. The radiative s-process timescale is quite 
long (a few 10* yr) compared with the convective s- 
process timescale (a few years). New stellar models 
of low mass stars (1 <M/M0 <3) were computed in 
the next 2 years, strengthening the new scenario of the 
radiative ^^C burning (Straniero et al. 1997). Gallino 
et al. (1998) carried out the full nucleosynthesis calcu- 
lation based on these new models. The mass and the 
profile of ^^C within the pocket were treated as free 
parameters. The Gallino ST case (corresponding to 
4 X 10^^ Mq of ^^C) provided the best reproduction of 
the main s-process component in the solar system. In 
particular, the Rb/Sr and ^^Zr/^°Zr problems as well 
as the energetic feedback problem found by Bazan and 
Lattanzio were removed. It is worth to note that the 
new stellar scenario for the s process was, in a certain 
sense, much more complex than previously believed. 
In fact, the different layers within the pocket undergo 
different neutron exposures, because of their different 
temperature and ^^C abundance. These differences 
are averaged when the s-rich pocket is engulfed into 
the convective zone generated by the subsequent ther- 
mal pulse, but due to the partial overlap with previ- 
ous convective zones, the freshly synthesized material 
is also mixed with the material processed by previous 
s-process episodes. In addition, the marginal activa- 
tion of the ^■^Ne(a, n)^®Mg may furtherly modify the 
heavy element distribution. As pointed out by Arlan- 
dini et al. (1999), the classical s-process analysis based 
on the simple exponential distribution of the neutron 
exposures is clearly unable to describe such a more 
complex nucleosynthesis scenario. 

More recently, Busso et al. (2001) and Abia et al. 
(2002) showed how the ratios between the flrst (Is), 
the second (hs) and the third (Pb) s-proccss peaks are 
sensitive to the mass of ^"^C in the pocket. This effect 
is illustrated in Figure 1. On this basis, Gallino and his 
team concluded that a certain spread in the ^"'C mass 
is required to explain the observed spread of [hs/ls] 
and [Pb/hs]. 

6 Conclusions 

In the last decades several hypotheses have been ad- 
vanced to explain the origin of the ^"'C pocket and 
quantitatively predict its mass and chemical profile. 
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and the surface overabundance of heavy elements with 
A>80 are significantly larger. 




[Fe/H] 

Figure 1: Variation of the spectroscopic index 
[hs/ls] as a function of [Fe/H] for different mass 
of ^'^C in the pocket (see text for the definition of 
the ST case). 



The most promising one was that of Hcrwig et al. 
(1997, sec also Hcrwig 2000), where a parameterized 
convective overshoot was adopted, based on prescrip- 
tions derived by 2D hydrodynamical calculations of 
stellar convection. They obtained ^^C pockets with 
mass of the order of 2-4x10"'' M© that is 10 to 20 
times smaller than the ST case defined by Gallino et 
aL (1998). Then, Langer et al. (1999) investigated the 
possibility of rotational induced mixing, while Denis- 
senkov & Tout (2003) the effect of a weak turbulence 
induced by gravity waves. 

Let us discuss in more detail the problem of the 
physical description of the convective boundary at the 
time of the third dredge-up. When the convective en- 
velope penetrates the H-exhaustod core, a steep vari- 
ation of the composition takes place at the convective 
boundary. The mass fraction of H sharply drops from 
about 70%, within the convective envelope, to zero, in 
the underlying radiative layer. The discontinuity in 
the composition induces a sharp variation of the ra- 
diative opacity and, in turn, an abrupt change of the 
radiative gradient. In these conditions, the precise lo- 
cation of the convective border (i.e. the limit of the 
region fully mixed by convection) becomes highly un- 
certain. A small perturbation causing a further mixing 
is amplified on a dynamical timescale, the radiative 
gradient in the radiative stable zone rises up and the 
convective instability moves toward the interior. This 
situation is commonly encountered in stellar model 
computations at the time of the second and the third 
dredge-ups (Becker and Iben 1979, Castellani, Chieffi 
and Straniero 1990, Frost and Lattanzio 1996, Mar- 
coni, CastcUaru and Straniero 1998). If the effect of 
such an instability may be maarginal in the case of 
the second dredge-up, this is not true for the third 
dredge-up. Two are the main consequences: i) deeper 
dredge-ups and ii) earlier occurrence of the first TDU 
episode. In other word, the minimim core mass for 
the occurrence of the TDU is lower. The observable 
consequences are many and important. The C star lu- 
minosity function is shifted toward lower luminosities 
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Figure 2: The boundary of the convective enve- 
lope during the third dredge-up. Lower panel: the 
exponential decline of the convective velocity and 
the sharp pressure gradient. Upper panel: chem- 
ical composition in the transition region between 
the convective envelope and the radiative core. 



A more realistic description of the convective bound- 
ary than that usually adopted in extant stellar evo- 
lution code is, in principle, required. Instead of a 
well defined spherical surface, as obtained when the 
bare Schwarzschild's criterion is used, the transition 
between the full-radiative core (i.e. unmixed) and the 
fuU-convective (i.e. fully mixed) envelope most likely 
occur in an extended zone where only a partial mixing 
takes place, so that a smooth and stable H profile may 
form. The evaluation of the actual extension of this 
transition zone and the degree of mixing there occur- 
ring would require sophisticated and reliable hydro- 
dynamical tools. At present, however, the inclusion 
of such kind of algorithms in the hydrostatic stellar 
evolution codes is probably a mission impossible. So 
far, only limited, even if enlightening, hydrodynami- 
cal investigations have been carried out (see e.g. Her- 
wig et al. 2006). Nonetheless, a different approach 
may be followed. When the convective envelope pen- 
etrates the H-cxhausted core, the average convective 
velocity at the inner border (ii, which is proportional 
to the difference Vrod — Vad) rises up (see Figure 2). 
Note that this occurrence confirms the instability aris- 
ing at the convective boundary during dredge-up. In 
order to control such an instability, in our most recent 
calculations of low mass AGB stars (Straniero et al. 
2006, Cristallo et al. 2009), wo have assumed that at 
the convective border the average convective velocity 
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drops smoothly to 0. In particular, we impose an ex- 
ponential decline of v, whose steepness is given by a 
free parameter to be calibrated in some way: 

« = «.ceexp(^-^) , (1) 

where d is the distance from the formal convective 
boundary (as defined by the Schwarzschild's criterion), 
vtce is the average velocity at the most internal con- 
vective mesh, Hp is the pressure scale height and /3 is 
a free parameter. A similar approach has been already 
followed by Herwig (1997), who assumed an exponen- 
tial decline of the diffusion coefficient at the convective 
border. 

The effects of the use of this algorithm on the phys- 
ical and chemical properties of the dredge-up are illus- 
trated in Figure 2. Note that in most cases since vtce is 
nearly zero at a convective border, the resulting tran- 
sition zone is extremely thin. On the contrary, during 
a dredge- up episode, owing to the growth of Vbce, the 
transition region, where v declines smoothly to zero, 
becomes quite extended. 




0.579 0.58 0.581 0.579 0.58 0.581 
M/Mj, M/Mg 



Figure 3: The formation of the C pocket. The 
sequence of panels A-B-C-D shows the evolution of 
the chemical composition in the transition zone be- 
tween the H-rich envelope and H-exhausted core. 
The various Unes represent: H (crosses), ^^C (dot- 
ted), (solid) and i^'N (dashed). Panel A: the 
third dredge up is just occurred and the convec- 
tive envelope is receeding; panel B: the tempera- 
ture and the density increase, the synthesis of ^'^C 
starts; panel C: where the number of protons is 
larger, some ^''N is also produced; panel D: the 
-'^^C and the "'^^N pockets are fully developed. 

As a result, a stable H profile is left in the He- 
rich intershell region after the third dredge-up, whose 
sharpness may be changed by acting on the v decline 
parameter /3. It is in this zone that, during the inter- 
pulse, the ^^C pocket forms. The total mass of 



depends on the /3 value. On the basis of numerical 
experiments made by changing the v decline param- 
eter we found i) the total mass of material dredged 
up increases with /3 and ii) there exists a maximum 
mass of '^^C in the pocket, as obtained for (3 about 
0.1, roughly corresponding to the ST case found by 
Gallino et al 1998 (Figure 3). For larger value of /3 too 
much protons are mixed into the transition zone, so 
that ^"N, rather than ^^C, forms. 




Figure 4: The surface composition at the AGB 
tip of 3 stellar models with solar metallicity and 
M = 1.5, 2.0 and 3.0 Mq, blue, black and red 
line, respectively, as obtained by means of our lat- 
est version of the stellar evolution code, which in- 
cludes a full network fro H to Bi. 



An important remark concerns the mixing algo- 
rithm. Our tests show that when a difi^usive scheme is 
adopted, as in many stellar evolution codes, the max- 
imum resulting ^^C mass is too small to allow a suffi- 
cient s-process nucleosynthesis in the pocket (see Her- 
wig 2000). On the other hand, it has been recently 
argued by Meakin & Arnett (2007) that the mixing 
occurring through the convective boundaries is not a 
diffusive process. As a matter of fact, '^^C pockets as 
large as those claimed by Gallino et al (1998) can be 
obtained if the degree of mixing between two (mesh) 
points depends linearly on the inverse of their recipro- 
cal distance, as in our time dependent mixing scheme 
(see Straniero et al. 2006). 

We have calculated new models of low mass stars 
by coupling the FRANEC stellar evolution code with 
a full network, from H to Bi. In such a way, we have 
been able to self-consistently derive the evolution and 
the nucleosynthesis along the whole AGB. The result- 
ing compositions at the AGB tip for 3 models of so- 
lar metallicity are shown in Figure 4. The first 
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pockets, those formed after the first 2 or 3 thermal 
pulse episodes, are quite large (in mass). Later on, fol- 
lowing the general shrinkage of the He-rich intershell, 
also the mass of the ^'^C pocket decreases. So, the 
resulting heavy element distribution at the surface is 
mainly determined by the first few ^^C pockets. Note 
that, owing to the overlap of the various He convec- 
tive zones, the imprint of these first bigger ^'^C pock- 
ets is conserved up to the last dredge-up episode. At 
high metallicity and in stars with initial mass smaller 
than 3 Mq, the temperature in the He convective zone 
never attains values sufficiently large for the substan- 
tial activation of a second neutron burst driven by the 
^■^Ne(a, n)^^Mg. Decreasing the metallicity, thermal 
pulses become stronger and generate larger temper- 
atures, so that at Z — 0.0001 the additional contri- 
bution of the ■^^Ne(a, n)^^Mg to the neutron capture 
nucleosynthesis is rather important (Cristallo et al. 
2009). 

Let us conclude by noting that the theoretical sce- 
nario built up so far, after about fifty years of intense 
investigations, appears to be firmly supported by the 
observation of abundances in AGB stars undergoing 
the third dredge-up and belonging to different stellar 
populations. 
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